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ABSTRACT 

We present a study of the X-ray properties of the nearby face-on Scd spiral galaxy MlOl 
based on recent XMM-Newton observations. In this third and final paper in the present series, 
we focus on the spatial and spectral properties of the residual emission, after excluding bright 
X-ray sources with Lx > 10'^^ erg s^^. Within a central region of radius 10' (21 kpc), the 
X-ray emission broadly traces the pattern of the spiral arms, establishing a strong link with 
recent star formation, but it also exhibits a radial scale length of « 2.6' (5.4 kpc) consistent 
with optical data. We estimate the soft X-ray luminosity within the central 5' (10.5 kpc) 
region to be Lx w 2.1 x lO'^^ erg (0.5-2 keV), the bulk of which appears to originate 
as diffuse emission. We find a two-temperature thermal model best fits the spectral data with 
derived temperatures of 0.20lo pi keVand0.68l[5:^^ keV which are very typical of the diffuse 
components seen in other normal and starburst galaxies. More detailed investigation of the X- 
ray morphology reveals a strong correlation with images recorded in the far-UV through to 
V band, with the best match being with the U band. We interpret these results in terms of a 
clumpy thin-disk component which traces the spiral arms of MlOl plus an extended lower- 
halo component with large filling factor. 



Key words: 

galaxies: individual: MlOl 



1 INTRODUCTION 

The hot (T ~ 10® K), X-ray emitting plasma observed in the disk 
of our Galaxy is an important factor when considering the energy 
balance and enrichment of the interstellar medium (ISM). Ener- 
gised by the winds of young massive stars and supemovae, its 
properties relate to the supernova rate, the evolution of supernova 
remnants (SNRs) and the impact of young star clusters on their 
environment. In normal spiral galaxies like the Milky Way, hot 
gas produced by multiple supernova explosions may escape from 
the d i sk via galactic fountains and/or chimney s dshapiro & Field 
1 19761 : lBregmatJ[l980l : |n orman & Ikeuchi|[T989h resulting in a dy- 
namic corona that is bound to the galaxy. This gas then cools, con- 
dens es and falls back to the disk, possibly as high- velocity clouds 
tereg man 1980). In contrast, the extremely high star-formation 
rates in starburst galaxies give rise to hotter (T > 10^ K) large- 
scale metal-rich outfl ows from the starburst regions in the form 
of "superwinds" (e.g.. lHeckman. Armus. & Milevlll990l) . Such su- 
perwinds are energetic enough to fully escape the influence of the 
galaxy and are potentially import ant for the enrichment and heating 
of the intergalactic medium (e.ig.. lMartin. Kobulnickv & HeckmanI 
l2003l) . 

Unfortunately, the opacity of the galactic plane to soft X- 
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rays greatly complicates the study of the extent, physical properties 
and filling factor of the X-ray emitting plasma in our own galaxy. 
However, observations of nearby spiral galaxies with configura- 
tions ranging from face-on to edge-on can circumvent many of the 
problems inherent in studies of the Milky Way. The Einstein obser- 
vatory provided the first detection of ultra-hot out-flows from the 
disks of edge-on starburst galaxies ( .Watson et al.. 1984 ). although 
only upper limits were obtained on t he diffuse emission produced 
in the disks of no rmal spiral galaxies dSregman & Glassgoldlll983 : 
iMcCammon & S anders 1984). The improved throughput and spa- 
tial resolution of subsequent ROSAT PSPC/HRI observations re- 
vealed u nambiguous evidence for truly di ffuse emission in spiral 
galaxies jCui et al.iri996l : lRead et al.lll997h . More recently, diffuse 
structures have been observed with the XMM-Newton and Chan- 
dra observatories in s amples of both edge- on (Pietsch et a l. 2001 
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2004b) and face-on spiral galaxies 1 Kuntz et al.1 



2003; Tvle r et alj|2004) . In the face-on systems, diffuse X-rays 
were found to correlate with the nuclear regions, trace recent star 
formation within the spiral arms and correlate spatially with Ha 
and mid-infrared emission jTvler et aL I l2004h . In the survey of 
edge-on systems by Strickland et al. (2004a,b), extra-planar dif- 
fuse emission was detected in all starburst galaxies and one normal 
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Table 1. Details of the thi-ee XMM-Newton observations of Ml 01 . 



# 


Observation ID 


Date 


Filter 


Pointing 


co-ordinates 


Useful exposure (ks) 






(yyyy-mm-dd) 




RA (J2000) 


Dec (J2000) 


pn 


MOS 


1 


0104260101 


2002-06-04 


Medium 




+54° 20' 24" 


25.7 


36.7 


2 


0164560701 


2004-07-23 


Medium 


14''03'"32.3'' 


+54° 21' 03" 


26.2 


28.9 


3 


0212480201 


2005-01-08 


Thin 


14''03'"32.7'' 


+54° 21' 02" 


27.7 


25.8 



spiral galaxy, the extent of which correlated with estimates of the 
rate of star formation within the disks. 

Clearly, if the purpose is to image the X-ray emission over a 
wide bandwidth, one should preferentially select galaxies in direc- 
tions where the line-of-sight Galactic hydrogen column density is 
low. In this context, the Scd supergian t spiral MlOl (NGC 5457) , 
for which iVj^ = 1.1 x 10^" cm"^ jPickev & LockmanllT990h . 
is an ideal candidate for a face-on study. A diffuse X-ray compo- 
nent was first hinted at in the central reg ions of MlOl with Ein- 
stein at energi es in the range 0.2-1.5 keV jMcCammon & Sanders! 
1 1984 : Trinchi erietalJll99(]h. S ubsequent ROSAT PSPC observa- 
tions l lSnowden & Pietschf 1995 ) revealed extensive diffuse X-ray 
emission in the 0.1-1 keV band, corresponding to thermal emission 
in the temperature range 1-3 x lO'' K. The dominant soft compo- 
nent was spatially peaked towards the centre of the galaxy but could 
also be traced out to a radius of at least 7' (15 kpc at a distance of 
7.2 Mpc). It was argued that this soft emission could not originate 
solely in the disk as this would require a filling factor greater than 
unity. After correcting for absorption intrinsic to MlOl, the total 
diffuse luminosity was estimated to be at least 10*° erg s~^. 

In the ROSAT HRI study of IWang et al"] ( Il999h . low-surface- 
brightness X-ray emission was detected both in the central region 
of MlOl and in the vicinity of bright HII regions in the spiral 
arms of the galaxy. The improved spatial resolution of the HRI 
(~ 5" ) allowed the detection and subtraction of point sources 
down to a luminosity of 4 x lO^'^ erg s^^, after which the total 
0.5-2 keV luminosity of the residual emission within a radius of 
12' was estimated to be 9 x 10 ^^ ergs~^. These authors agreed 
with ISnowden & PietscH ( 1 19951) that faint discrete X-ray sources 
were likely to account for only a fraction of the observed flux, im- 
plying that the bulk of the soft emission is truly diffuse in nature. 
The X-ray spectra of the emission in the central region from both 
the ROSAT PSPC and ASCA GIS were fitted simultaneously, yield- 
ing a soft thermal component with a temperature of kT ~ 0.2 keV, 
plus a hard power-law continuum most plausibly explained as the 
integrated emission of the population of hard X-ray binary sources 
(XRBs) in MlOl. 

Further prog ress was made in the Chandra study of 
iKuntz et alj feOOsI) . using data from the 8' square ACIS-S3 chip 
from the first ~100ks Chandra observation of MlOl. The superb 
spatial resolution of Chandra allowed the removal of point sources 
down to luminosities of ~ 10'^® erg s~^, leaving residual emission 
of luminosity 2.3 x 10^^ erg s"^ in the 0.45-1.0 keV band, 16% of 
which was estimated to come from unresolved XRBs and dwarf 
stars. The bulk of the soft emission was found to trace the spiral 
arms and was also spatially correlated with Ha and far ultra-violet 
(UV) emission, physically linking the soft X-ray emission to re- 
gions of on-going star formation. The X-ray spectrum of the diffuse 
emission was well-described by a two-temperature thermal plasma 
with fcr=0.2/0.75 keV, with a large covering factor implying that a 
significant fraction of the softest component originates in the halo 
of MlOl. 



This is the third in a series of papers repo rting the re- 
sults o f XMM-Newton observations of MlOl. In Ijenkins et al.l 
( l2004al) (hereafter Paper I), the spectral and timing pr operties of the 
brightest X-ray sources were investigated, whereas in'lenkins et alj 
( 2005b ) (hereafter Paper II) we reported on the full catalogues of X- 
ray source detected in MlOl by XMM-Newton. In the current paper 
we focus on the morphology and spectral properties of the "diffuse" 
X-ray component in this galaxy. The structure of this paper is as 
follows. First we describe the XMM-Newton observations and the 
methods used to construct images and spectra of the residual emis- 
sion after the exclusion of bright discrete sources (section[2l(. Next 
we investigate the likely composition, spectral properties and spa- 
tial morphology of this residual emission (section O. We then go 
on to discuss the implications of our results (section|4l and, finally, 
provide a brief summary of our conclusions (section lSTl. Throughout 
this p aper we assume a distance to MlOl of 7.2 Mpc Tstetson et alH 
1 19981) . implying that an angular scale of 1' corresponds to a linear 
extent of 2.09 kpc in MlOl. 



2 OBSERVATIONS & DATA REDUCTION 

Three XMM-Newton EPIC observations have been made of MlOl 
as summarised in Table [T] Our earlier analysis of the point source 
population in MlOl, presented in Papers I and II of this se- 
ries, utilised only the first observation performed in 2002. How- 
ever, for this paper, we have taken advantage of two additional 
observations carried out by XMM-Newton in 2004/5 as a Tar- 
get of Opportunity (TOO) programme to follow up the outburst 
of t he ultraluminous supersoft X-ray source J140332.3 -I-542103, 
(see lKong, Pi Stefano, & YuarJl2004lMukai et al.ll2005l) . The first 
XMM-Newton observation was targeted at the nucleus of the 
galaxy with the result that the 30' diameter field-of-view of the 
EPIC cameras e ncompassed the full 28.8' (D25) extent of MlOl 
Vaucouleurs et al. 19_91). In the second and third observations 
the source J140332. 3-1-542103 was placed on axis, resulting in a 
3.3' eastward offset of the field-of-view with respect to the first 
observation. 

All datasets were screened for periods of high background by 
accumulating full-field 10-15 keV light curves. MOS data were ex- 
cluded during periods when the 10-15 keV count rate exceeded 0.3 
ct s~^, and the pn data were similarly screened when the pn count 
rate exceeded 3 ct (except in observation 1, which was more 
heavily contaminated by flaring, where the pn threshold was set at 
0.9 ct s^^). For the subsequent image and spectral analysis, sin- 
gle and double pixel events were selected for the pn (pattern 0-4), 
whereas single to quadruple events (pattern 0-12) were utilised in 
the case of the MOS data. 



2.1 Image processing and point source subtraction 

To study the morphology of the diffuse component, pn and MOS 
images and associated exposure maps were created for each ob- 
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servation in three energy bands: soft (0.3-1 keV), medium (1- 
2keV) and hard (2-6 keV). In each case the pixel size was set at 
4.35" X 4.35". The pn, MOS 1 and MOS 2 images in each of the 
three bands were first flat-fielded by subtracting a constant (non- 
vignetted) background particle rate (estimated from the count rates 
recorded in the corners of the CCD detectors not exposed to the 
sky) and then dividing by the appropriate exposure map. In the 
same process bad pixels, hot columns and spurious data along chip 
gaps were excised. The resulting data from the three observations 
were then mosaiced into flat-fielded pn and MOS (MOS I plus 
MOS 2) images for each band. Since the medium and hard band 
images showed no evidence for diffuse emission, here we focus 
primarily on the images obtained in the soft (0.3-1 keV) band. Fig- 
ure [Ha) shows an adaptively- smoothed soft-band image, obtained 
from the combined pn and MOS datasets, encompassing the cen- 
tral 20' diameter region centred on the optical nucleus of MIOI 
(RA 14'^ 3™ 12.55", Dec +54° 20' 56.5" ) . The amplitude scal- 
ing is logarithmic and adjusted so as to highlight the extended low 
surface brightness emission in MlOl. 

To facilitate the investigation of the contribution of bright 
sources to the total X-ray emission in MIOI, a master sourcelist, 
applicable to the three observations, was created. To this end, each 
of the new observations was analysed in turn using the source de- 
tection procedure described in Paper II (applied to the soft, medium 
and hard band data). Then, by comparing the resulting source lists, 
the 108 sources detected in observation 1 (Paper II) were supple- 
mented by 21 new sources from observation 2 and an additional 9 
new sources from observation 3, to give a combined source cat- 
alogue of 138 sources (see the Appendix for details of the new 
sources). 

The next step involved the construction of a spatial mask for 
the suppression of the bright X-ray sources in the central field of 
MIOI. Using the pn-l-MOS flat-fielded soft-band image, a soft-band 
source count rate, net of the local background, was determined at 
each source position; in practice, we extracted the count rate within 
a cell of 16" radius centred on the source position, subtracted the 
local background, and lastly applied a correction for the extent of 
the point spread function (PSF) beyond the source cell. We then 
selected the sources within II' radius of the nucleus of MIOI with 
count rates above 1.5 pn-FMOS ct ks~^ as a "bright source sam- 
ple" to be excised from the the soft-band image. Here our approach 
was to construct a model image in which each source in our bright 
source list was represented by a PSF sub-image centred on the 
source position and an intensity scaling commensurate with the 
measured source count rate. Then by applying a cut at a surface 
brightness threshold of 0.01 pn-l-MOS ct ks^^ pixel"^, we were 
able to construct a bright source mask - see Fig.fHb). A total of 65 
bright soft-band sources were removed within a radius of 1 1' of the 
nucleus of MIOI (25 sources within 5' ). 

Figure[TJc) shows the pn-l-MOS soft-band image of M 1 1 after 
applying the bright source mask with the colour table adjusted to 
highlight the distribution of the residual emissiorQ. 

The residual component is most prominent within the the cen- 
tral 5' radius region, although it can be traced out to 10' in the radial 



profile distribution shown in Fig.|2] There is a clear correlation with 
the spiral structure of the galaxy, as is evident from a comparison 
with the GALEX near-UV (NUV; 2310 A) band image of MIOI 
teianchi et aI.ll2005l : |PoDescu et alj2003) shown in Fig.fTJd). 



2.2 Spectral extraction of the diffuse component 

The images in Fig. [T] demonstrate the existence of an extensive 
residual component in the soft X-ray emission of MIOI. In order 
to investigate the spectral properties of this component we have 
extracted spectral data for the 5' radius region, centred on the po- 
sition of the nucleus. A total of nine individual diffuse emission 
spectra were derived from the appropriate event files (for the three 
EPIC instruments, over the three observations). In each case, re- 
gions around bright point sources were excluded using the mask 
described earlier (in practice this involved the use of circular exclu- 
sion cells of radius 12" — 53", as an approximation to the source 
mask used in the image analysis). 

Given the extended nature of the diffuse emission, the back- 
ground subtraction is more complicated than is usual with EPIC 
data. The use of other XMM-Newton data sets deriving from "blank- 
field" observations is, in our view, not a good option due to the 
variation from field to field in the soft X-ray background from the 
sky - clearly an important consideration when attempting to mea- 
sure a soft diffuse signal. Instead local background (LBG) spectra 
were extracted from each event file for the annular region from 5' to 
10' radius - see Fig. [TJa) - again with the regions around bright 
sources excluded. We also extracted particle background (PBG) 
spectra from the out-of-field-of-view corner regions in each EPIC 
camera. The true background beneath the residual MIOI emission 
in the central source region will be made up primarily of (vignetted) 
photons from the X-ray sky and (non-vignetted) particles, which 
comprise the instrumental background. The LBG spectra will also 
have photon and particle components, but, as the LBG region is 
further off-axis than the source region, these components are not 
in the correct ratio for accurate background subtraction (since the 
signal reduction due to the vignetting is larger in the LBG region). 
However, a very good model of the true background can be made 
from a combination of the LBG and PBG spectra, provided the 
contributions of the two are appropriately scalecQ to allow for the 
vignetting-corrected and actual areas of the central, LBG and PBG 
extraction regions. Using this approach, an appropriate composite 
background was determined for each central region spectrum to- 
gether with the corresponding ARE and RMF files. Finally for each 
dataset, the spectral channels were binned together to give a mini- 
mum of 20 counts per bin. 

Given the radial profile in Fig. |2] the above amounts to a dif- 
ferential measurement of the galaxy spectrum, which in effect com- 
pares the integrated signal in the central 5' region with the equiva- 
lent signal measured in the 5' - 10' annular region. 



^ Hereafter in this paper we use the word "residual" to describe the X-ray 
emission of MIOI after masking out the bright discrete X-ray sources. This 
residual component will include the truly diffuse emission, the integrated 
emission of discrete sources with count rates below our source-exclusion 
threshold and also some low-level contamination from bright sources due 
to the spillover beyond the masked regions of their extended PSFs - see 
Fig.Elb). 



^ The scale factors applied to the LBG and PBG spectra were respectively 
(Ac/Albg) X (Vc/Vlbg) and (A^/Apbg) x (1 - Vc/Vlbg)- Here 
Ac, Ai^BG and ApBG refer to the actual areas on the detector of the cen- 
tral, LBG and PBG extraction regions. Similarly Vc and Vx^bg are the 
average vignetting factors within the central and LBG regions. 
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Figure 1. (a) XMM-Newton soft-band (0.3-1 keV) image of the central 20' of MlOl. This flat-fielded image is constructed using pn and MOS data from the 
three XMM-Newton observations. The amplitude scaling is logarithmic with a dynamic range (maximum to minimum in the colour table) of 50. The position of 
the optical nucleus of MlOl is marked with a cross, (b) A "model" image constructed from a set of count-rate scaled PSF sub-images. A source mask, defined 
by the black contours, was produced by applying a cut at a surface brightness threshold of 0.01 pn-l-MOS ct ks^^ pixel^^. (c) The same X-ray image as 
shown in panel (a) with the regions confused by bright sources masked out. The amplitude scaling is logarithmic with a dynamic range of 10. (d) The GALEX 
near- UV(NUV:2310A) band image of MlOl on the same spatial scale as the X-ray data. The scaling is logarithmic with a dynamic range of 10. In all four 
panels the two large circles have radii of 5' and 10'. Figures submitted to arXiv as jpgs . 




3 RESULTS 

3.1 Luminosity and composition of the residual emission 

After the exclusion of the bright source regions and subtracting the 
mean sky surface brightness measured in the 5' -lO'annulus, the 
count rate in the central 5 region is 0.28 pn-^MOS ct s"^ (0.3-1 
keV). However, we need to apply three correction factors to obtain 
the "true" central 5' count rate. The first accounts for the loss of 
"diffuse" signal due to the source mask which blocks 12% of the 
central area. The second factor allows for the differential nature of 
the measurement {i.e., the fact that we extract the background from 
the 5' -10' annulus where the galaxy is still bright); on the basis of 
the radial profile determined earlier we estimate that this effect sup- 
presses the central count rate by ~ 19%. Finally, the spillage of the 
bright source signal beyond the masked region (due to the extended 
PSF) contributes about 9% of the measured residual count rate (as 
determined from the source model image shown in Fig.fTJb)). The 
correction of all three effects translates to a 28% upward scaling of 
the residual count rate. If we adopt the best-fitting two-temperature 
model described in section [T2l the count rate to flux conversion 
factor is 1 pn-l-MOS ct s"^ = 0.97 x 10"^^ erg cm"^ s"^ for the 
0.5-2.0 keV band (fluxes corrected for Galactic absorption). The 
equivalent factors for the 0.3-1 keV and 0.45-1 keV bands are 1.16 
and 0.91 x 10~^^ erg cm~^ (ct s"^)"'^ respectively. 

The implied luminosity of the central emission of MlOl after 
excluding bright sources is 2.1 x 10"^^ erg in the 0.5-2 keV 
band. However, we must take into account the fact that our source 
exclusion threshold is relatively high (roughly a fact or 10 higher 
than that employed in the equivalent Chandra study - iKuntz et al.l 
12003 ':). Chandra has shown that the log N - log S relation for rela- 
tively faint sources in the central region of MlOl (after correction 



for background AG Ns) is quite flat w ith a slope of —0.80 ± 0.05 in 
the integral counts jPence et al.l200lh . implying that bright sources 
provide the dominant contribution to the discrete source luminosity 
of the galaxy. Using the count rate to flux conversion factor quoted 
earlier, our source exclusion threshold of 1.5 pn-l-MOS ct ks^^ 
(0.3-1 keV) corresponds to a source luminosity of 9 x 10^® erg s~ ^ 
(0.5-2 keV), although this could be an underestimate for the harder 
spectral forms more typical of bright X-ray binaries. In fact, a 
source exclusion threshold of 1 x 10''^ erg s^^ appears to be a 
reasonable estimate based on the average count rate to flux con- 
version ob tained when we co mpare our measurements with those 
quoted in Pence et al] ( 12001 ) for the sub-set of common sources. 
As a check, at the latter luminosity threshold the Chandra source 
counts predict « 25 sources within a central 5' region, which is 
very consistent with the number of sources actually excluded in 
our study. Using the Chandra source counts, we find the integrated 
luminosity in discrete sources from our exclusion cut-off down to a 
factor 10 deeper (in effect the Chandra limit) amounts to about 20% 
of the residual signal measured by XMM-Newton. In other words, 
our estimate of the "diffuse" luminosity of MlOl (after allowing 
for the contribution of sources brighter than ~ 10^'' ergs~^) is 
1.7 X 10^^ erg (0.5-2 k eV) in good agreement with the value 
reported bv lKuntz et alj ( l2003i) based on Chandra measurements of 
the same region. 

3.2 Spectral properties of the residual emission 

The nine epic spectra were fitted simultaneously using a variety of 
spectral models, with identical model parameters applied to each 
of the data sets. Visual inspection of the spectra indicated almost 
no emission above « 1.5 keV. No acceptable fit could be obtained 



© 2005 RAS, MNRAS 000,[T]l9] 



An XMM-Newton view of M 101 - III 5 



using just a single hot plasma component. However, a good fit (re- 
duced x^=0.87 for 2986 de grees of freedom) was obtained with 
a model incorporating two thermal (solar-abundance) Mekal com- 
ponents, with absorption attributable to the hydrogen column den- 
sity of our Galaxy in the direction of MlOl (1.16x 10^" cm~^; 
iDickev & LockmarJ [1990). The derived temperatures for the two 
components were 0.20+o oi keV and 0.6Sto,ol keV (90% confi- 
dence limits assuming 1 interesting parameter). The data plus best 
fitting two-temperature model are shown in Fig [3] (In this figure 
we show the averages of the three pn and the six MOS spectra for 
clarity of display; however, the spectral fitting results quoted here 
pertain to the fits to the nine individual spectra). The absorption- 
corrected flux in the 0.3-1 keV band derived from the spectral fits 
was 3.0 X 10^^'^ erg cm^^ s^^ of which the hotter component con- 
tributes 18%. The corresponding values for the 0.5-2.0 keV band 
are 2.5 x 10"^"^ erg cm~^ s^^ with the hotter component contribut- 
ing 27%. When we allowed the metallicity of the two plasmas to 
vary from unity (i.e., solar values) we obtained a best fit of 0.044 
(0.038-0.055) solar for the lower-temperature component and 0.17 
(0.10-0.50) solar for the hotter component (with a reduced for 
the fit of 0.81). However, this determination of strongly sub-solar 
abundances is most probably an artifact of our attempt to fit rel- 
atively low spectral resolution data pertaining to a complex multi- 
phase plasma with a simplistic spectral model, a pr oblem which has 
in fact been well documented in the literature (e.g., Strickland et al. 
L2OOO; Wang et al. 2001). As noted bv .Kuntz et a l. (2003), the abun- 
dances in the hot plasma in the central region of MlOl are more 
likely to be near to solar rather than grossly sub-solar. 

We know that the measured XMM-Newton spectrum includes 
some residual contamination by the wings of the bright sources (es- 
timated to be 9% of the 0.3-1 keV count rate) and also includes 
a contribution from the underlying unresolved source population 
(estimated to contribute 20% of the 0.5-2 keV flux). Since the 
sources may have somewhat harder spectra than the diffuse emis- 
sion, presumably they contribute preferentially to the higher tem- 
perature spectral component. However, our derived temperatures of 
« 0.2 and ~ 0.7 keV are in very good agreement with the results 
of iKuntz et al. Illooi based on Chandra observations, suggesting 
that some fraction of the harder emission does have a diffuse ori- 
gin. In the following analysis we assume 50% of the measured hard 
emission is attributable to diffuse emission. We also scale the com- 
ponent normalisations obtained from the spectral fitting so as to 
correct for both the area of the source mask and the differential 
nature of the spectral measurement (as detailed in §3.1). 

One can infer mean physical properties of the hot diffuse gas 
once some assumptions have been made regarding the geometry 
of the diffuse emission. Here we assume that the gas is contained 
within a cylindrical region of radius of 10.5 kpc in the plane of the 
galaxy (matching a 5' angular extent) with a half- width perpendicu- 
lar to the plane of 0.5 kpc (representing the extent of a putative shal- 
low halo). Using this volume, the derived emission measure rjniV 
(where r] is the 'filling factor' - the fraction of the total volume V 
which is occupied by the emitting gas) can be used to infer the mean 
electron density n^. For the 0.2 keV plasma we find Tie ~ 0.003 
rt~^^'^ cmT^ compared to a value of ~ 0.001 vf^^"^ cvoT^ for the 
0.7 keV plasma, implying that these two components are in rough 
pressure balance. The thermal energy residing in each of the plasma 
components is comparable and totals Eth ~ 6 x 10^^ r;^'^^ erg. The 
cooling is dominated by the line emission with radiative cooling 
timescales of t^^^^ ~ 1.8 x 10 77 ' yr and t^ard ~ l-S'? Gyr 
for the cooler and hotter components respectively. 



3.3 Morphology of the residual emission 

As previously noted, the spatial distribution of the residual soft 
X-ray emission in MIOI shows a good correlation with the opti- 
calAJV emission of the galaxy, at least over the central 5' region. 
We have investigated this correlation using the FUV (1530 A) 
and NUV (2310 A) im ages of MlOl from the GALEX pipeline 
jMorrissev et al.l 1 20051) , together with images recorded by the 
XMM-Newton Optical Monitor (OM) in the UVWl (f« 2800 A), 
U, B and V filters (Mason etal. 2001). The latter were derived 
from data taken during observation 1 (UVWl) and observation 3 
(U,B,V) and subject to the standard SAS pipeline processing. 

To aid the comparison, the soft X-ray image of Fig.[TJc) was 
compressed in both coordinates by a factor 4 (to give 17.4" x 
17.4" pixels) and then lightly smoothed. The UV and optical im- 
ages were first rebinned into the same image format as the original 
soft X-ray data, blurred by a mask representative of the on-axis 
XMM-Newton soft X-ray PSF and then further compressed in an 
identical way to the soft X-ray data. 

Fig. m shows a detailed comparison of the 0.3-1 keV surface 
brightness measured in the central 5' radius region of MIOI with 
the corresponding NUV, U and V band images. A striking similar- 
ity between the soft X-ray and U band images is immediately obvi- 
ous. In particular, the "S-feature" which delineates the inner spiral 
structure of MlOl is very prominent in both cases. The fall-off in 
the surface brightness towards the edges of the 5' field is also very 
similar in both wavebands, with the north-south asymmetry evident 
at the level of the lowest soft X-ray contour level also present in the 
U band image. Interestingly, the correspondence becomes slightly 
less strong when the soft X-ray image is compared to both the NUV 
and V band images. The NUV image is clearly totally dominated 
by bright HII regions and young star associations which together 
delineate the spiral arms of the galaxy, whereas in the soft X-ray 
image the general rise in surface brightness towards the nucleus 
has more relative impact. It would seem highly unlikely that these 
morphological trends could be induced by the increase in extinction 
as we move from the U band into the NUV. In contrast, comparison 
with the V band data veers in the other direction with the soft X-ray 
image showing more structure outside the central 2' zone than the 
optical image. 

The above trends are further illustrated in Fig. [5] which 
shows the pixel-by-pixel correlation between the soft X-ray sur- 
face brightness and the optical/UV surface brightness in the vari- 
ous wavebands, (here we use unsmoothed data with a pixel size of 
34.8" X 34.8" , i.e. a factor of 8 compression of the original im- 
ages). The correlation, as measured by the Pearson product moment 
correlation coefficient, R, is at a maximum in the U band. 



4 DISCUSSION 

After excluding bright point sources (with Lx > W^'' erg s"^) 
we measure the residual soft X-ray luminosity of the inner (10.5 
kpc radius) disk of MIOI to be Lx ~ 2.1 x 10^'' erg s"^ (0.5- 
2 keV). Up to 20% of this flux may originate in relatively bright 
point X-ray sources (with Lx down to 10^^ erg s~^), which re- 
main unresolved in the XMM-Newton observations. A further 5% 
of the measured luminosity may be attributable to intermediate lu- 
minosity sources (with Lx down to 10"^"* erg s~^) and, at very faint 
levels, the integrated e mission of dwarf st ars may account for an- 
other 10% contribution dKuntz et al.ll2003h . Nevertheless, it seems 
that the bulk of observed X-ray emission is likely to be diffuse in 
origin. 
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Figure 3. The EPIC spectra of the residual X-ray emission component in MlOl after the exclusion of bright sources, (a) The pn spectrum averaged over 
the three observations, (h) The combined MOS spectrum averaged over the three observations. In both cases the solid line corresponds to the best-fitting 
2-component thermal spectrum with the dashed lines showing the individual contributions of the 0.2-keV and 0.7-keV components - see text. The residuals 
with respect to the best-fitting model are also shown. 



Figure 4. (a) The pn+MOS soft-band (0.3-1 keV) image of the central 5' radius region of MlOl. The scaling is linear The contour levels increase in steps of 
0.3 pn-l-MOS ct ks-i pixel"! ^^^^ ^ starting level of 0.9 pn+MOS ct ks^^ pixel"!, ffoj GALEXNUW image of the same region overlaid with the soft 
X-ray contours. The scaling is linear in count pixel" ^ units. ( c) The XMM-Newton OM U image of the same region overlaid with the soft X-ray contours. The 
scaling is linear in count pixel"! j-^j jj^g XMM-Newton OMV image of the same region overlaid with the soft X-ray contours. The scaling is linear in 

count pixel"! units. Figures sTj]3].nitted to arXiv as jpgs. 
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The EPIC spectrum of the residual emission can be modelled 
as the sum of two thermal plasma components at temperatures of 
~ 0.2 keV and « 0.7 keV with the harder component contribut- 
ing 27% of the flux in the 0.5-2 keV band (18% in the 0.3-1 
keV band). These spectral characteristics appear to be quite typi- 
cal of th e diffuse co mponents seen in norma l and starburst galax- 
ies (e. g., iFratemalTet al. 2002; Jen kins et al. ||2004b; Jenkins et al., 
l2005ah and, albeit with some caveats, consiste nt with spectral mea- 
surem ents for MlOl reported in the literature dSnowden & PietschI 
1 1995. : Wang 1999 ; Kuntz et alj|2003h . 



Earlier studies have demonstrated that for normal spiral galax- 
ies, the tota l soft X-ray lum inosity c orrelates with the total infrared 
luminosity teead & Ponm an 2001 ; S trickland et alj2004bl) . the lat- 
ter serving as a direct indicator of the current star formation rate. 
For the diffuse X-ray components, this is in line with the mecha- 
nism whereby the mechanical energy input from supernovae and 
the winds from young massive stars serves to heat the galactic ISM 
to million-degree temperatures. In the case of MlOl the required 
rate of energy i nput is readily supplied by two supernovae per cen- 
tury jMatonick & Fesen 1997 ). if the efficiency of the conversion 



of the mechnical energy of the supernova explosion into X-rays is 
10^^ /E %, where E is the energy of the supernova in erg. 

The linkage of X-ray emission to recent star formation is fur- 
ther strengthened by the finding that in a number of galaxies the 
X-ray morpholo gy in the spiral arm s matches that seen in the mid- 
infrared or Ha ( lTvIeretal.ll2004 h. A recent Chandra study has 
demonstrated that, in MIOI, the diffuse X-ray emission traces the 
spiral arms and is c orrelated with both the Ha and far-UV emission 
jKuntz et alj2003l) . In the same study there was no evidence for any 
significant trend of X-ray spectral hardness with radiufl implying 
that the 0.2- and 0.7-keV spectral components are closely linked 
with one another and are both associated with star formation. 

In the present paper we have investigated this correlation in 



^ We have subsequently checked the XMM-Newton data for spectral trends 
by extracting images in the 0.3-0.8 keV and 0.8-1.2 keV sub-bands. The ra- 
dial profiles of the residual X-ray emission in these sub-bands are consistent 
with the results obtained for the broader 0.3-1.0 keV band (see Fig.|2) and 
within the en'ors, exhibit a constant ratio for radii in the range I'-IO'. The 
XMM-Newton observations thus confirm that there is no significant radial 
dependence of spectral hardness in this Galaxy. 
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more detail by contrasting the X-ray morphology with that ob- 
served in far-UV through to V band. We find that the correlation 
with the X-ray morphology is quite strong for all bands consid- 
ered (FUV, NUV, UVWl, U, B and V), although the best match is 
obtained with the U band data. One interpretation of this result is 
that, in broad terms, there are two underlying components present 
in the X-ray surface brightness distribution, namely a smooth rise 
towards the centre of the galaxy as reflected in the radial profile of 
the emission (Fig. |2]( plus a superimposed X-ray spiral arm struc- 
ture. In moving from the FUV through to the V band, clearly the 
balance shifts between the extreme population I components (dis- 
tributed in the spiral arms) and a more intermediate population (as- 
sociated with the galactic disk) and it is in the U band that this 
balance happens to best match the X-ray morphology. 

In the disks of spiral galaxies star-formation is thought to be 
triggered by the passage of a spiral density wave through the ISM. 
The formation of massive young stars results in a large UV/far-UV 
flux with reprocessing in the immediate environment subsequently 
giving rise to associated Ha and mid-IR emission. All of these 
serve as tracers of the spiral arms of the galaxy. Population synthe- 
sis models predict that following a star-formation burst, the peak in 
the UV/Hof production will last ~ 3 x 10® yr, whilst the most mas- 
sive stars complete their life cycle. The rate of energy deposition 
into the ISM from the subsequent supemova e rises at this time and 
remains fairly constant for ~ 3 x 10^ y r jLeitherer & HeckmanI 
ll995l : ICervino. Mas-Hesse & Kunthll2002h . Given the lag between 
the peak in UV/Ha emission and the X-ray heating one might pre- 
dict a spatial offset between the diffuse X-ray emission and the 
other spiral tracers, althoug h as yet there is n o clear observational 
evidence for such an effect jTvler et al]|2004b . 

In the case of M 101, we assume a flat rotation curve be- 
yond a few kpc of the galaxy centre a nd a rotational velocity 
Vrot ~ 200 km s"^ teosma, GosT^Al len 1981). We further 
assume that the spiral pattern corotates with the disk material at 
3 radial scale lengths, i.e. at r « 15 kpc. Then at r « 7.5 kpc a 
delay of 3 x lO'^ yr translates to a rotational lag of 24° or about 
3 kpc in a direction perpendicular to the spiral arm. The effective 
pixel size in our multi-waveband correlation is 35" or 1.2 kpc, im- 
plying that a delay in the X-ray production following the passage 
of a spiral density wave might just be measurable in the current 
data. Unfortunately the restricted signal to noise and limited clarity 
of the spiral arm features in the X-ray image mitigates against the 
detection of such an effect. Earlier we estimated a radiative cool- 
ing timescale for the more prominent 0.2-keV plasma component 
to be t^gff ~ 1.8 x 10* yr for a filling factor 77 ~ 1. Clearly this 
timescale is inconsistent with the presence of reasonably narrow 
spiral features in the soft X-ray image. A very small filling factor, 
e.g. T) ~ 0.01, is required to match the properties of X-ray spiral 
arms in MlOl, suggesting a very clumpy thin-disk distribution. In 
fact a reasonable scenario is that the spiral arms in X-rays are de- 
lineated by a combination of truly diffuse emission, possibly in the 
form of hot gas bubbles and superbubbles, plus contributions from 
individual S NRs and concentrati ons of unresolved discrete sources. 
As noted by iKuntz et^ ( l2003h . the available constraints suggest 
that superbubbles in the disk of MlO l may hav e similar properties 
to the Galactic Loop I superbubble jEgger & A schenbach 1995), 
although they are not individually resolvable even in Chandra data. 
Where there is no confinement by chimneys or similar structures 
in the ISM, energy losses arising from adiabatic expansion of the 
hot gas in the disk into the lower halo of MlOl may also help lo- 
calise the spiral arm component {nb. in MlOl the sound speed of 



the s oft component is c omparable to the escape velocity from the 
disk - iKuntz et alj|2003h . 

The second spatial component of the X-ray emission con- 
sidered above broadly follows the distribution of the optical light 
attributable to the intermediate disk stellar population. The inte- 
grated X-ray emission of dwarf stars will presumably follow a disk 
distribution but this is unlikely to represent a substantial contri- 
bution (10% overall in the 0.5-2 keV band but, given the typi- 
cal kT ~ 0.8 keV temperature of dwarf-star corona e, falling to 
f» 5% in the softer 0.3-1 keV band: lKuntz et alj2003h . Other con- 
tributions to this smoother distribution might be expected from su- 
pernovae in the interarm regions associated with disk population 
stars and X-ray emitting gas which has accumulated in a shal- 
low halo in MlOl as a result of galactic chimney/fountain activ- 
ity. In the case of the latter, gas which has cooled substantially 
due to adiabatic expansion may, as a result of a frozen-in non- 
ionization equilibrium, exhibit an emission spectrum dominated by 
oxygen lines, which mimics that ex pected from a 0.2 keV plasma 
jSreitschwerdt & Schmutzleij|l999t) . The fact that the X-ray sur- 
face brightness increases towards the centre of MlOl is most natu- 
rally explained in terms of an increase in activity per unit disk area, 
with the filling factor of the extended z-height component possi- 
bly approaching unity near the centre of the galaxy. With the latter 
assumption, the implied electron density and pressure in the inner 
halo of MlOl is comparable to that inferred for the centre of Galac- 
tic Loop I and the Local Hot Bubble. 



5 CONCLUSIONS 

This is the third and final paper in a series presenting the results of 
XMM-Newton observations of the nearby face-on Scd spiral galaxy 
MlOl. Here we focus on the spatial and spectral properties of the 
galaxy, when bright X-ray sources with Lx > 10^^ ergs~^ are 
removed using an appropriate point source mask. 

The residual soft X-ray luminosity of the central (10.5 kpc ra- 
dius) region of MlOl was measured as Lx ~ 2.1 x 10^^ erg 
(0.5-2 keV), the bulk of which appears to originate as diffuse emis- 
sion. We find a two-temperature model best fits the spectral data 
with the derived temperatures of 0.20±0.01 keV and 0.68lo 34 
keV, typical of the diffuse components seen in other normal and 
starburst galaxies. 

In line with earlier studies, we find that the observed X- 
ray surface brightness distribution is well correlated with images 
recorded in optical/UV wavebands. In particular the detection of 
spiral arms in X-rays establishes a close link between the X-ray 
emission and recent starformation. Closer investigation suggests 
that the X-ray morphology may comprise both a spiral arm com- 
ponent and a smoother disk component. In spiral galaxies, star 
formation is thought to be triggered by the passage of a spiral 
density wave through the ISM. In principle one might observe a 
lag between more immediate star-formation indicators such as the 
UV/far-UV flux or Ha and the heating of the interstellar medium to 
X-ray temperatures as result of subsequent supernova. The XMM- 
Newton data show no evidence for such an effect but it may be 
observable in future high resolution and high signal-to-noise data. 
Whereas the spiral arm component, on the basis of the radiative 
cooling timescale of the 0.2-keV plasma, may be deduced to have 
a clumpy, low-z distribution, the smoother disk component may 
represent longer-lived X-ray gas with a relatively large filling fac- 
tor, which has been transported to the lower halo of MlOl via the 
galactic chimney/fountain mechanism. 
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Future progress in understanding the X-ray properties of nor- 
mal spiral galaxies will no doubt follow from very intensive studies 
of individual galaxies (e.g., the recent 1 Ms observation of M 101 
with Chandra - Kuntz et al. in preparation) and by applying appro- 
priately optimised analysis procedures to samples of nearby galax- 
ies, as is now possible using archival datasets. 
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APPENDIX: THE POINT SOURCE CATALOGUE 

Within the D25 circle of M 101, 108 sources were detected in obser- 
vation 1, 91 in observation 2 and the same number in observation 3. 
Ninety sources were detected in at least two of the observations. 
For completeness, Table |2] lists the additional 21 sources detected 
in Observation 2 and another 9 sources detected in Observation 3. 
The format of the table is identical to that of Table 1 of Paper II and 
the numbering scheme is a continuation of the original list of 108 
sources. 

The most notable new source to appear in the second and 
third observations is source #122, wh ich corresponds to the super- 
soft transient P98 jMukai et al]|2003h . The results from the XMM- 
Newton and Chan dra observations of thi s object are reported in 
iKong et al.l ( l2004t) and lMukai et"ai] ( |2005|) . 
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Table 2. Additional sources detected within the D25 circle of MlOl in observations 2 (#109-129) and 3 (#130-138). 



Src 


XMMU 




pn count rate (count ks 


-1) 


MOS cou 


nt rate (count ks ^ } 


F V 


Lv 


HRl 


HR2 






(") 


S 


M 


H 


s 


M 


H 










(1) 


(2) 


(3) 




(4) 






(5) 




(6) 


(7) 


(8) 


(9) 


109 


J140213.8+542158 


1.95 


3.1±0.9 


0.1 ±0.3 


1.7±0.8 


0.8±0.2 


0.9±0.2 


0.7±0.2 


1.58±0.27 


0.98±0.17 


-0.47±0.12 


0.21±0.16 


110 


J140237.7+542733 


2.98 








0.4±0.2 


0.9±0.2 


0.3±0.2 


0.90±0.30 


0.56±0.19 


0.41±0.24 


-0.52±0.28 


111 


J140238.0+541215 


2.31 








0.5±0.2 


1.0±0.3 


0.3±0.2 


1.11±0.32 


0.69±0.20 


0.34±0.21 


-0.49±0.26 


112 


J140246.7+541454 


2.01 


0.0±0.3 


2.8it0.6 


1.3 ±0.5 


0.0±0.1 


0.7±0.2 


0.4±0.2 


1.00±0.21 


0.62±0.13 


1.00±0.14 


-0.32±0.15 


113 


J140251.4+540901 


2.93 








1.2±0.3 


1.2±0.3 


0.4±0.2 


1.56±0.39 


0.97±0.24 


0.00±0.18 


-0.49±0.26 


114 


J140257.9+542655 


2.07 


l.OibO.4 


1.4it0.5 


1.5±0.5 


0.2±0.1 


0.5±0.2 


0.8±0.2 


1.27±0.22 


0.79±0.13 


0.25 ±0.20 


0.16±ai5 


115 


J140302.0+541336 


2.00 


1.8ib0.5 


0.8it0.4 


0.0±0.2 


0.6±0.2 


0.2±0.1 


0. 1 ±0.1 


0.40±0.12 


0.25±0.08 


-0.43±0.17 


-0.67±0.41 


116 


J140304.0+541040 


2.07 


2.8ib0.8 


1.2it0.5 


0.0±04 


0.5±0.2 


0.8±0.2 


0.1 ±0.2 


0.64±0.18 


0.39±0.11 


-0.07±0.15 


-0.76±0.28 


117 


J140306.9+541007 


2.28 








0.3 ±0.2 


0.0±0.1 


1.2±0.3 


1.79±0.47 


1.11 ±0.29 


-1.00±0.89 


1.00±0.22 


118 


J140308.3+543229 


2.01 








1.4±0.3 


0.0±0.1 


0.0±0.0 


0.60±ai7 


0.38 ±0.10 


-1.00±ai3 


0.00±0.00 


119 


J140320.0+542034 


2.93 








0.7±0.2 


0.1±0.1 


0.0±0.1 


0.35±0.14 


0.22±0.09 


-0.70±0.24 


-1.00±1.36 


120 


] 140324.0+542336 


2.58 


2.S±0.5 


04±0.3 


0.0±0.1 








0.34±0.11 


0.21 ±0.07 


-0.75±0.17 


-0.83±0.69 


121 


J140327.5+541909 


1.82 


2.1±0.5 


1.1 ±0.3 


0.8±0.3 


0.4±0.1 


O.S±0.1 


04±0.1 


0.88±0.14 


0.55±0.09 


-0.14±0.14 


-0.12±0.16 


122 


J140332.4+542103 


1.73 


S.5±0.7 


0.1±0.2 


0.6±0.3 


0.9±0.2 


0.0±0.0 


0.0±0.1 


0.55±0.10 


0.34±0.06 


-0.98±0.04 


0.64±0.45 


123 


J140334.3+540930 


1.95 


4.2±0.9 


2.7±0.7 


1.9±0.8 


1.0±0.3 


0.5±0.2 


0.9±0.3 


1.91±0.32 


1.19±0.20 


-0.26±0.13 


0.04±0.17 


124 


J140354.1+541 108 


2.20 








1.4±0.3 


0.0±0.1 


0.1 ±0.2 


0.79±0.27 


0.49±0.17 


-1.00±0.17 


1.00±1.74 


125 


J140355. 8+542100 


1.68 


6.3±0.8 


2.2±0.S 


0.7±0.4 


0.9±0.2 


0.7±0.2 


0.1±0.1 


1.03±0.14 


0.64±0.09 


-0.38±0.08 


-0.62±0.16 


126 


J140357.2+ 54101 1 


2.11 


S.Sil.l 


2.3±0.7 


1.1±0.7 


0.6±0.3 


04±0.3 


1.0±0.3 


1.67±0.32 


1.03 ±0.20 


-0.38±0.14 


0.02±0.21 


127 


J140401.2+542344 


2.59 


2.3±0.5 


0.0±0.2 


0.3±0.4 








0.44±0.23 


0.27±0.14 


-1.00±0.16 


1 .00± 1 .09 


128 


J140429. 8+542235 


2.19 








0.3±0.2 


0.8±0.2 


0.6±0.2 


1.25±0.30 


0.78±0.19 


0.43±0.23 


-0.18±0.22 


129 


J140435.4+542101 


1.59 


12.9±1.3 


0.0±0.2 


0.0±0.2 








1.41±0.20 


0.88±0.12 


-1.00±0.03 


-1.00±1.00 


130 


J140229.9+542243 


1.54 


14.0±1.4 


0.0±0.3 


0.0±04 








1.44±0.26 


0.89±0.16 


-1.00±0.04 


0.00±0.00 


131 


J140230.0+541611 


2.05 


3.6±1.1 


4.4±1.0 


3.2±1.1 


1.2±0.3 


1.2±0.4 


0.7±04 


2.30±0.42 


1.43±0.26 


0.07±0.14 


-0.20±0.17 


132 


J140243.6+542007 


2.38 


2.1±0.6 


0.4±0.3 


04±0.4 








0.48±0.24 


0.30±0.15 


-0.70±0.25 


0.04±0.71 


133 


J140250.7+542857 


3.01 


2.5±0.7 


0.0±0.1 


1.3±0.7 








1.00±0.37 


0.62±0.23 


-1.00±0.09 


1.00±0.17 


134 


J140300.6+540958 


2.24 


4.3±1.0 


2.0±0.7 


1.0±0.8 








1.29±0.45 


0.80±0.28 


-0.36±0.18 


-0.36±0.38 


135 


J140312.2+541754 


2.19 


3.8±0.7 


0.3±0.3 


0.1 ±0.2 








0.50±0.14 


0.31 ±0.09 


-0.84±0.13 


-0.45±0.76 


136 


J140314.4+542132 


2.79 








0.6±0.2 


0.2±0.1 


0.2±0.1 


0.57±0.19 


0.35±0.12 


-0.49±0.25 


-0.15±0.48 


137 


J140355.8+540856 


2.10 


4.0±0.9 


2.9±0.8 


1.3±0.7 








I.57±0.44 


0.97±0.27 


-0.17±0.I7 


-0.39±0.27 


138 


J140415.6+540948 


1.75 


7.2±1.1 


4.5±0.9 


1.5±0.9 


2.0±0.4 


1.0±0.3 


0.6±0.3 


2.12±0.31 


1.31±0.19 


-0.26±0.09 


-041±0.17 



(1) source number; (2) XMMU source designation (J2000 coordinates); (3) la error radius (including a 1.5 arcsecond systematic error); 
(4 & 5) source count rates in soft (0.3-1 keV), medium (1-2 keV) & hard (2-6 keV) bands for the pn and MOS cameras, with the 
significant source detections (> 4(t) highlighted in bold; (6) source flux in units of 10^^'* erg cm^^ s^^ in the broad (0.3-6 keV) band; 
(7) source luminosity in units of 10''** erg in the 0.3-6 keV band (assuming a distance to MlOl of 7.2 Mpc); (8 & 9) soft (HRl) and 
hard (HR2) hardness ratios. 
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